Control of Towing Kites for Seagoing Vessels 

Michael Erhard and Hans Strauch 



Abstract — In this paper we present the basic features of 
the flight control of the SkySails towing kite system. After 
introduction of coordinate definitions and basic system dynamics 
we introduce our model used for controller design and justify 
its main dynamics with results from system identification based 
on numerous sea trials. We then present the controller design 
which we successfully use for operational flights for several years. 
Finally we explain the generation of dynamical flight patterns. 

Index Terms — Aerospace control, Attitude control, Feedfor- 
ward systems, Wind energy 



I. Introduction 

THE SkySails system is a towing kite system wliich allows 
modern cargo ships to use the wind as source of power 
in order to save fuel and therefore to save costs and reduce 
emissions |[T|. The SkySails company has been founded in 
2001 and as main business offers wind propulsion systems for 
ships. Starting the development with kites of 6-10 m-^ size the 
latest product generation with a nominal size of 320 m^ can 
replace up to 2MW of the main engine's propulsion power 
Among the marine application of kites there is a strongly 
increasing activity in using automatically controlled kites fS), 
lia, f4l, f5l, f6l, fJ], f8l and rigid wings f9l, fW\ in order 
to generate power from high-altitude wind |11J. Since 2011 
the company's second business segment SkySails Power also 
develops and markets systems for generating power from high- 
altitude wind. Therefore the design of control systems for kites 
in these geometries has become a growing field of theoretical 
O, Uni, ini, IISJ and experimental ll9l, 

Il20l . ED research efforts. 

The main components of the SkySails system are shown 
in Fig. [T] The core of the propulsion system is the towing 
kite steered by the control pod situated under the kite. The 
towing force is transmitted to the ship by a high-strength 
synthetic fiber rope. Additionally a launch and recovery system 
is installed aboard the ship One key component of the 
flight control system is the main steering actuator in the control 
pod applying deflections to some kite lines leading to curve 
flight. A control system consisting of distributed computers 
preprocesses data from various sensors at a rate of 10 Hz 
and performs the flight control algorithm which calculates the 
steering command applied to the main actuator. An integrated 
graphical user interface allows for operation of the system 
by the crew whereas for research and development purposes 
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Fig. 1. The BBC SkySails with a 320 m^ towing kite. 



prototyping and testing toolchains can be connected via special 
interfaces. 

The paper is organized as follows: First we introduce the 
basic system and coordinate definitions. We then focus on 
the main dynamics and develop a model specially suited for 
controller design. After justification of the main law of the 
model with experimental data we present our controller design 
discussing design considerations and controller performance 
measurements. We complete the article with the explanation 
of pattern generation. 

n. Basic System and Coordinates 

In this section we give a mathematical description of the 
considered system. It is worth mentioning that we deal with 
a constrained system which shows a completely different 
dynamics compared to free flying parafoils |22|. The basic 
system is sketched in Fig. |2] Compared to the real system 
we make use of some simplifications. The flexible rope is 
substituted by a rigid rod which also is parallel to the kite 
yaw axis Cyaw and all masses are neglected. At first view this 
seems to be an unrealistic simplification, but the usual mode 
of operation is the highly dynamical pattern flight leading 
to line forces big compared to system masses. Therefore 
inertia effects or free flight situations, where the towing line 
is no longer stretched out due to gusts or wave induced ship 
motions, are infrequent. Although consideration of these issues 
becomes important at a certain point when bringing the system 
to higher perfection, a detailed description of the solution to 
these off-nominal situations would go beyond the scope of this 
paper 

We would like to start with a demonstrative and introductive 
example to motivate the idea of the chosen coordinate system 
ip, d and to explain the basic dynamics. Imagine a kite flying 
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Fig. 2. Definition of coordinates for tlie considered system. Tlie right- 
lianded coordinate system is defined by the basis vectors ex,ey, with Sx 
in wind direction and pointing downwards with respect to gravity. The kite 
position is parameterized by introducing the spherical coordinates </3 and i? 
(for a more precise definition see 1(2)). The kite axes are labeled as roll ejoji, 
pitch Epitch ™d yaw eyaw ■ This corresponds to the definition usually used in 
aerospace applications with roll axis parallel to forward and yaw axis parallel 
to down directions respectively. Note that the yaw vector eyaw is defined by 
the position of the kite assuming it is constrained to the origin by a rigid rod. 
Thus orientation of the kite is represented by the single angle tj}. Detailed 
vector definitions are given in the appendix of the paper. 
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Fig. 3. Flying a kite in a fictitious wind tunnel on a space station would 
be instructive in understanding the basic equations. While the neutral flight 
is stationary at any arbitrary position like for example A, a constant ip 
due to steering will lead to a circular orbit B. The diameter is a function of 
the tp magnitude. See also Fig. [4] for further illustration. 



in a wind tunnel experiment conducted on a space station 
which means the absence of gravity. We further assume the 
free manoeuvrability of the kite unrestricted by obstacles 
like ground or ceiling. From this mental picture it would be 
natural to chose a coordinate system with symmetry axis in 
wind direction. A kite with its roll axis 6^011 antiparallel to 
the wind direction would stay at a some arbitrary, stationary 
position like e.g. 'A', see Fig. |3] However, once a deflection 
is commanded, the kite will go 'down' and orbit around 
permanently on a circular path 'B'. The diameter of the steady 
state circle is a function of our coordinate which in turn is 
determined the by commanded deflections. The smaller the 
diameter of the circle the faster flies the kite and thus the 
more force will be generated. The corresponding dynamical 
equations will be discussed in detail in section [Till 

We would like to emphasize that the above sketched model 
is our controller design model while we use for various 
other development and test issues a sophisticated simulator 
model including multi-body dynamics which also captures the 
aerodynamical effects more comprehensively (comparable in 



structure to ll23l ). Yet, we suggest above model based on this 
imaginary space station experiment because it is specifically 
suited for the controller design purpose. It neglects gravity 
effects and by this a new symmetry axis is gained. The main 
benefit of this symmetry compared to other coordinate systems 
lH], [24 1, 1 25 1 is the resulting simplicity of equations of motion 
which allows to approach the feedback and guidance design 
task to a large extent analytically or semi-analytically at least. 
Further on it provides a straight forward way of describing 
flight patterns. 

Our choice of the design model led to a controller structure 
to be presented in section [V] This controller structures turned 
out to work quite effectively in numerous sea trials. 

We close this section by definitions for the quantities used 
in the following. The coordinate system is shown in figure [2] 
For a constant line length L the state of the kite is defined by 
the three angles ip, d and V'- With respect to the basis vectors 



we obtain for the kite position x: 

COS1} 
L I sin If sin {} 

— cos (f sin {} 



(1) 



The kite axes are denoted as ej-on (roll or longitudinal), 
Gpitch (pitch) and Cyaw (yaw). An explicit definition of these 
vectors is given in appendix [A] 

For a description using rotation matrices one would start 
with a kite at position Le^ with roll-axis in negative z- 
direction ej-on = —Sz and then apply the following rotations: 
—tfj about X, I? about y and finally ip about x. This transfor- 
mation reads; 



R = R,{p)Ryi^)R,i-^b)■ 



(2) 



One could interpret the angle ?/; as orientation of the kite 
longitudinal axis with reference to the wind. For a given 
kite position x (parameterized by (p and z?) the reference 
orientation ip — corresponds to the minimum of the scalar 
product (eioihetr) obtained when turning the kite fixated at 
this position x around its yaw axis eyaw A nonzero value 
represents a kite orientation obtained by a rotation of about 
the yaw axis Cyaw starting at this reference orientation. 

III. System Dynamics used for Design 

For verification and other development purposes we use 
a full dynamics simulation containing a multi-body model, 
an aerodynamic database and parameters adopted to results 
of sea trials. In this section we would like to present the 
main relations and the dynamics of a complementary model 
specifically tailored for the design of the controller The 
detailed derivation steps are summarized in the appendix. The 
equations of motion for d and (p read: 

, tani?\ 
cos tp — — (3) 



Va_ 

L 



E 



(4) 



L sin "d 

Thus the dynamic is mainly controlled by the angle -tp- Further 
quantities are the air path speed Va, the towing line length 
L and the glide ratio E. As already pointed out and to 
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Fig. 4. Zenith positions for neutral fliglit witli low force. With ^/i = the 
parameter can be freely chosen determining the positions (drawn as white 
kites). For constant »/) the resulting i? converges as shown by the gray kites 
flying on the marked trajectory. 



be reasoned in the next section we neglect acceleration and 
gravity effects in order to obtain these simple equations of 
motion which allow for the following interpretations. 
First we get a steady state angle i9o for constant by 



1^0 ("0) — arctan(i?cos^/'). 



(5) 



A further result is the dependence of the (steady state) air path 
velocity Va on the value of and the ambient wind speed vq. 



Va = wq-E'cos'!? = -i;o-E'cos(arctan(i? 008-0)) 



(6) 



which is the key issue for pattern generation as we can use i/; as 
a tuning knob to control Va and the force which is proportional 
to v'^ accordingly. It is worth mentioning that for practicing 
the sport of kite surfing, the content of (j6|l is crucial in order 
to control forces: kite surfers know the angle as position in 
the so called 'wind window' and the deeper they fly their kite 
into this 'wind window' — i.e. they decrease this angle — the 
more traction force they get and vice versa. 

A special case is the so called zenith position for i/; = 0, 
a neutral flight situation which generates only low forces and 
thus is used mainly for launching and recovering the kite. 
From (j5| we get i9q — arctan(i?) and ip can be used as free 
parameter to determine the neutral flight position (compare 
Fig.|4j|. 

While all equations up to here follow straightforward from 
kinematic reasoning, albeit the choice of the coordinate system 
was not 'typical', the dynamic response of the kite to a steering 
deflection S is claimed to be 



(7) 



where g is the proportionality factor We would like to draw 
special attention to this turn rate law (|7| and will show in 
the following that it can be justified by measured data to a 
surprising high degree. Therefore it is a key issue for the 
cascaded controller approach where it constitutes the dynamics 
of the inner loop. 

Finally we would like to point out that due to its motion 
on a spherical surface an inertial sensor measures a turn rate 
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Fig. 5. Flight trajectory under computer control for the bang-bang flight 
experiment. Angles were measured by tow point sensors on the ship deter- 
mining the direction of the towing rope. Note that azimuth and elevation 
values correspond to a different coordinate system which is not defined in 
this paper, but still show up the figure-eight shape clearly. 



tpm about the yaw axis Cyaw different from the derivation 
Tp = dijj/dt. The rotation measured by the pod sensors can 
be calculated by transforming the dynamics represented by R 
into the pod coordinate system by applying Comparing 
the rotation operation fix with R-R^ yields: 



■0™ = ^ 'f cos -d. 



(8) 



By consideration of flight trajectories one can assure oneself 
that the second term of this equation usually is small compared 
to the first and thus can be neglected in the first instance and 
treated as a correction to the controller design later 

IV. Justification of the Turn Rate Law 

In this section a justification of (|7]i is given based on 
numerous experiments showing the strong proportionality. It is 
worth noting that (j7| is confirmed by measurements to a high 
degree even in disturbed sea trial conditions. The key issue of 
these experiments is to perform bang-bang flights which will 
be presented and discussed in the following. 

The excitation of the system for the identification is per- 
formed in the following way: We apply a constant steering 
command +5o to the system. The system will respond with a 
positive yaw rate (0™ > 0). When reaching a certain threshold 
0m > 00 the corresponding opposite steering deflection —6a 
is commanded leading to a decrease of ?/;,„. Falling below 
the negative threshold 0„i < ^jq, the primary deflection Sq is 
commanded again. The schedule of the bang-bang experiments 
is as follows: the human pilot flies the kite into a high zenith 
position with tp^iQ and hands over the steering to the computer 
based control system which performs the described algorithm. 

A typical flight trajectory is shown in Fig. |5] The bang- 
bang steering leads to a figure-eight-pattern and for typical 
parameters the air path speed and thus the size of the pattern 
increase because the kite flies down to smaller elevation angles 
'0. The human pilot only has the task to supervise the flight 
and overtake manual control before the system runs into the 
danger of overload or bounces against the water surface. 
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Fig. 6. Measured data of a bang-bang flight. Yaw rate ipm as function of 
the air path speed multiplied by deflection VaS to justify jT}. The parameter 
g is obtained by the shown linear fit. 
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Fig. 7. Comparison of steering command 5 with yaw rate ipm and the yaw 
rate divided by the air path velocity 4>m/{g Va)- Note the increasing rates are 
due to increasing air path speed Va while going down into the wind window 
(compare Fig. [sj. 



In Fig. |6] data points for one experiment run are shown. 
This experiment was performed in 2006 using a 20 m^ kite. 
Repeating this experiment with different (5o leads to similar g 
values and thus proves the validity of (|7]i. 

Although the linear dependence can be clearly identified in 
Fig.|6]it is even more convincing to present the data in the time 
domain as shown in Fig. |7] Here we compare the time-series 
of the steering command with the turn rate of the kite. The 
trapezoidal shape is due to the finite steering velocity of the 
control pod. The resulting measured yaw rate i>m shows an 
increase which results from the increasing air path velocity 
Va over the experiment. The yaw rate divided by gva is 
also plotted in order to compare it to the steering command. 
Although a lot of perturbations affect these experiments we 
observe an excellent correlation. This analysis justifies the 
validity of (|7]i to a high degree and recommends its usage 
as a key role for the controller design. 

At this point we would like to classify and review these 
bang-bang experiments in a historical context. Following the 



textbook approach in classical system identification we did 
a lot of identification flights in the years 2005 and 2006 
using separate batch runs in order to characterize the steering 
behavior of our kites at the various operating points. The 
experiments were quite cumbersome: perturbations from wind 
gusts can be comparatively large compared to periodic exci- 
tations caused by the deflection commands and the air path 
velocities were difficult to tune for the different batch runs 
and desired operating points. As we had to evaluate data from 
different days with changing environmental conditions and 
drifting flight properties of our kites we solely could suspect 
the validity of law (|7|. But once we switched to a bang-bang 
flight strategy the real law shows up clearly. This is obvious 
as one bang-bang experiment varies the parameter Va over a 
whole range, in one flight alone, lasting only a few tens of 
seconds and therefore plays a trick on perturbations. 

We would like to conclude this section by giving a extended 
version of (|7]i which also takes into account the effect of the 
gravitational force on the turn rate and reads: 



tpni = gVaS + M 



COS On sin 'ipn 



(9) 



The quantity 9g denotes the angle between eyaw and the 
etc -Ley-plane and ipg the angle between eroii and the ex-Lcy- 
plane. 

Because a steering deflection could be regarded as a kite 
force component into pitch direction epitdi subsequently lead- 
ing to a yaw rate, the gravity force, projected onto the pitch 
axis by cosOgSintpg, should have the same effect. We have 
shown in this section that the yaw rate is proportional to Va S. 
This can be attributed to a side force proportional to S 
from aerodynamical and design considerations. Transferring 
this reasoning to the mass term, which is independent from 
Va, we expect a factor of M/va between the 'gravitational' 
side force and the yaw rate. The constant M includes system 
masses and kite characteristics. As 1\I is positive for our kites 
we get an instable behaviour and thus have to stabilize V'm by 
active control. 

As for the usual operation point of dynamic flight we have 
{gVa) '3> (M/va) and thus the second term of (j9|, which we 
call 'mass term', can be neglected for the design of the linear 
feedback law only to be introduced as correction term via a 
feedforward path to the controller structure. 

As our operational flights under autopilot conditions (dur- 
ing the dynamic flight modes) utilize similar bang-bang like 
commands, we can use the discussed identification scheme as 
a standard tool to establish or check the kite parameter during 
normal operation. During flights a recursive least-square algo- 
rithm ll26l runs in order to determine the system parameters 
g and M on-line. We monitor these values for changes which 
may indicate upcoming material failures in advance and can 
adopt controller parameters on-line, if necessary, to increase 
robustness. 

After we have convinced ourselves of the validity of (j7]i, 
we will now present the controller design which is strongly 
influenced in its structure by the discussed law. 
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Fig. 8. Cascaded controller approach for i/i control implementing the model 
following structure. Detailed diagrams for the blocks [FF^], [C^\, [FF^\ 
and [C^] are shown in figs. |13| |12| p| and llOl The Controller calculates a 
steeling command & from the input value ^/j^using a measured yaw angle 
ipm and yaw rate ipm- 



V. Controller Design 

Most of time the kite is operated in a highly dynamic regime 
where the air path speed can easily vary by a factor of up to 3- 
5 within some seconds and the deflection command can change 
by more than ±60%. The classical way to approach such a 
controller design task would be to use a controller structure 
which specifically aims at time varying, non-linear systems. 
Non-linear dynamic inversion or non-linear model predictive 
control could be such candidates. Gain scheduling based on 
linearized plant models along the trajectory is a further alterna- 
tive. Actually we tried the latter approach in the beginning but 
where not satisfied with the achievable robustness. The core 
of the problem is governed by the fact that it is difficult to 
execute a classical modelling approach as usually performed in 
aerospace application. Such an approach would be based on 
an aerodynamic database covering the full dynamic regime. 
Performing wind tunnel tests would be expensive and for our 
larger kites with 160 m^ area downright impossible. We also 
learned that a controller structure where the kite trajectory is 
given as a set-point directly into a unique controller block, 
which then directly computes the deflection command, lacks 
robustness due to the above mentioned modelling issue. 

Instead we settled for a separation of the overall controller 
task into 'guidance' and 'control'. Section [Vl| describes how 
the guidance algorithm computes 'tpg, which is then the input 
to the controller However the major distinguishing feature of 
our autopilot, compared to other approaches we found in the 
cited literature, is the cascaded controller structure which is 
based on the model following principle as shown in Fig. [8 

The basic idea is to reflect the separation of dynamics (7]) 
and kinematics, as given by t/j = dip/dt, adequately in the 
controller structure. We can see two cascaded loops in Fig. |8] 
The inner loop gets a commanded rate ipg as an input and 
computes the deflection command 5. The outer loop has V's 
as input and commands ips to the inner loop. 

From (j7]i we recognize that the dynamics from deflection 
to yaw rate can be viewed as a proportional plant, = 
S, where K^p — gva denotes the gain. Of course we 
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Fig. 9. Details of block [FF^\ of the cascaded controller (see Fig. jsj. 
From input value ips the feedforward value Sfj and are calculated by 
using a steering pod model consisting of a limiter and a ratelimiter Note the 
translation of rate to command and back via division and multiplication by 
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Fig. 10. Feedback block [C^] of the cascaded controller (see 
The controller mainly consists of a PI feedback on the yaw rate. Note tee 
division by , which is a function of the air path speed, thus introducing a 
nonlinearity into the feedback by transforming the rate command from the PI 
controller to the command portion Sf^jj which is routed to the steering pod. 



have to realize that is not constant but a function of 
the air path velocity. We take care of this by employing a 
feedforward/feedback structure which implements the model 
following principle: in the feedword term [FF^^] we compute, 
in an open loop fashion, the deflection command Sff necessary 
to achieve ips . The feedback control [C^] only acts when there 
is a remaining control error due to external disturbances or due 
to unmodeled plant dynamics. An appropriate delay block is 
necessary in order to capture all the delays from command to 
actual execution. 

With such a feedforward/feedback structure we can accom- 
modate the dependence of the gain from the air path 
velocity. Fig. |9] provides details of the feedforward controller 
block. In line with the idea of the model following principle 
the feedforward block is not limited to linear equations and can 
accommodate any system description. The extended version of 
the turn rate law (|9]i 



V'r, 



sin Tpg 



(10) 



can easily be considered. The principal idea is to invert it in 
order to compute the necessary deflection 5//. In addition the 
block also contains nonlinear elements, like limiters on angle 
and rate, in order to capture limited pod steering speed and 
other constraints. The mass term from (fTOll is introduced via 



M cos 6g sin i/jg 



(11) 



Fig. 10 provides the details of the feedback controller block: 
As the plant has a proportional character the feedback structure 
is of proportional/integrator (PI controller) type. The output of 
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Fig. 11. Experimental results for the inner loop. In the upper plot the 
measured yaw rate ipm, is compared to ips and i/'c- The lower plot compares 
feedforward Sjf and feedback ^fbi^ controller output signals for a typical 
dynamical flight situation. 
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Fig. 12. Feedback block [C^] of the cascaded controller (compare Fig.[8|. As 
the dynamics to be controlled is mainly an integrator, a proportional feedback 
has been chosen. 



the PI controller is divided by thus taking the velocity de- 
pendence into account. A lowpass is added for noise rejection 
of the measurement. 

By providing Fig. [TT] we illustrate how effective the inner 
loop actually works. We show how a couple of repeating 
flight patterns look at the inner loop level. The feedforward 
command 6ff moves between ±60% of the available deflection 
range as can be seen in the lower graph of the figure. The 
feedback command (5fbk hardly needs to correct control errors 
due to unmodeled plant dynamics. Actually this figure is just 
an alternative account to Fig. [6] in proving the good fit of the 
turn rate law. 

As seen from the outer loop the inner loop has dealt with 
the aerodynamically influenced part of the dynamics. It now 
remains for the outer loop controller to achieve a desired i/jg. 
The division into feedforward and feedback parts is kept. Figs. 
[T2] and [T3] provide further details. 

As the plant characteristic is of integrating nature a feed- 
back law with proportional character [C^], augmented by a 
lowpass, is sufficient. As before a limiter is also employed. 
The feedforward block [FF.ff,] 
features. 



(see Fig. 13 1 has more elaborate 



Note that the feedforward term features an internal 
feedback loop. The basic idea is to shape the commanded i/jg, 
even if it is a jump, in such a way that it corresponds to the 
actual response capabilities of the control pod and kite. We 
achieve this by employing an internal loop from t/jg to tp'^. 
Note that the two limiters are crucial in shaping the final tp'^ 



Fig. 13. Feedforward block [FF^] details of cascaded controller (see Fig.jsj. 
The steering pod model is included as combination of limiter and ratelimiter 
Note that even for step inputs on i/>s the algorithm computes feedforward 
ip ff and '0c in a way which is consistent with the capability of the overall 
dynamics. This is achieved by using an inner feedback loop embedding the 
steeling pod model. 
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Fig. 14. Flight results for the outer loop. The upper plot compares the 
response of the measured angle rpm to a given rectangular ips- The curve 
4>c is computed by the internal loop of the feedforward block [FF^] using a 
model of the steering pod (compare Fig. |13| . 



evolution. Furtherm ore w e have a nonlinear function inside 
the loop: f{x) = \J 26pX. We will stop short in deriving the 
details of this special feature. Instead we will illustrate it with 



Fig. 14 The upper part shows how the input command i/;^ is 
shaped into the command ipc by the non-linear feature of the 
internal feedback loop, tpc corresponds to a an actually flyable 
ipc pattern, tp ff shows the corresponding necessary rate pattern 
which is fed into the inner loop. 

As a conclusion we will summarize the basic design prin- 
ciples of our controller. The first feature is the separation of 
the dynamics of deflection to rate (j7| from the kinematic of 
rate to angle. This separation allows to introduce non-linear 
elements (mainly limiters) at the appropriate places. That way 
we achieve a shaping of our commanded signals such that 
they corresponds to the limitations of the complete chain from 
software command over control pod steering to kite movement. 

A second characteristical feature is the feedforward/feed- 
back separation which allows to decouple non-linear elements, 
as for example the mass term, from the feedback. The feed- 
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Fig. 15. Geometry for pattern generation. The figure-eight pattern is guided 
by two states (1) and (2). Transitions between those are triggered by the 
conditions Pi: ip <ipQ—ipa and P2: ip> ipo~ipa- The corresponding sequence 
is shown in Fig. [16] 



back loops can then be designed within the realm of linear 
control theory. Only proportional or integral dynamics remain 
which can be handled in a classical way. Due to the feed- 
forward/feedback separation the selection of the closed loop 
bandwidth of the two loops is more concerned with achieving 
sufficient stability margin than with achieving performance in 
terms of fast response because this is already mastered by the 
feedforward. 

The next chapter will illustrate the generation of the -ips 
command from the desired kite trajectory. Although it could 
be perceived as just a further cascaded loop we treat it more 
like the 'guidance' feature of classical aerospace applications. 

VI. Pattern Generation 

In this section we describe the dynamic flight mode which is 
used to generate traction force by flying dynamical patterns in 
order to obtain high air path speeds and forces. The algorithm 
utilizes the presented controller design by providing the input 
value ijjg- 

In order to explain the main principles we would like to 
review the space station experiment of sectionjll] In this model 
a constant value of -05 = +'i/'o or —ij^vi leads to a circular orbit 
clockwise or counterclockwise dependent on the sign of ijjs 
and the obtained force can be easily controlled by the value 
of \4>s\- For the purpose of line force generation in our space 
station experiment this would finish our design effort — but 
in our application we are not able to fly circular orbits as the 
kite would crash onto the water surface. Thus the solution is 
to turn around at certain points of the orbit and fly back and 
forth. 

The resulting trajectory of such a scheme is shown in 
Fig. [15] The underlying algorithm is similar to those of 
the bang-bang experiments in section IV A constant value 
■0s = +-00 is commanded until point 'Pi' is reached (at 
f < 'Po — fa) triggering the command tps = ^ 0o until point 
'P2', then triggering the former value t/jg = +ipo and so on. 
The timing of tps is depicted in Fig. 16 

It is acceptable to apply a square signal to as the 
controller design contains an internal model, which leads to 
a calculation of ip^ which is feasible to fly. We would like to 
mention that every curve involves a decrease of air path speed 
as -d increases and hence a decrease of traction force. As a 




Fig. 16. Pattern generation sequence toggling between the two states (1) 
and (2) when conditions Pi : ip < ipo — fa and P2 : > ipo — (fia are met 
respectively. The states directly result in a square signal to ^ which is shaped 
by the feedforward block [FF^] into •t(>c (compare sec. [Vj. The respective 
pattern geometry is displayed in Fig. |15| 



consequence curves should be commanded in an efficient way 
using appropriate deflections and manageable steering speeds. 
This is fulfilled by the controller design as it both utilizes the 
given curve deflection capability and the given design speed 
of the control pod for step transitions (see section [V|. 

As illustrated in Fig. 15 there are three parameters defining 
the trajectory. The parameter tpa determines the pattern size; 
the parameters ipQ and determine the center point of the 
pattern. The value ipQ can be freely chosen within a certain 
range by the operator or an overlying algorithm in order to 
optimize the force component pointing in forward direction of 
the vessel. The value -do can not be tuned directly but indirectly 
by the ipo value which is the tuning knob for the air path speed 
Va and hence force. Nevertheless (|6]l does not hold in its simple 
way and could by improved by using (-0) — J dt\tl;{t) \ instead 
of i/jQ in order to estimate i9 or resulting forces. 

As it is cumbersome to predict the exact wind situation at 
flight altitude in any case we make use of another approach for 
force control: we use an outer loop force controller evaluating 
the height of the force peaks while flying the figure-eight 
pattern which provides a feedback value for ipQ. Details of this 
control law as well as of the supervision mechanisms during 
operative flight and start procedure of the pattern are important 
and interesting issues each but would exceed the purpose of 



this paper In Fig. 17 we present the trajectory of some eights 
and show the corresponding time series of the angles ip„i, i/jg, 
ipc as well as of the air path and wind speeds. 

For sake of completeness we would like to note control 
strategies for the neutral flight mode introduced in section |ll] 
The determining equation is Q and the issue is to control ip 
to the given set value tpo. We mainly use a linear controller 
based on a classical PID controller in order to compute 4's 
from ipm~^o- 

VII. Concluding Remarks 

Dealing with a complex system in a demanding environment 
we could present several further topics on our control system 
which are closely related to the discussed contents. These 
topics have been omitted in the previous sections because 
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Fig. 17. Flight test results illustrating pattern generation using a 160m^ kite 
at an operational towing line length of 300 m. The trajectory for a typical 
figure-eight pattern is shown in the upper plot (compare Fig. |15( . The two 
plots in the middle show corresponding curves for ipm, 4>c and 4'm (see 
Fig. |16| . The lower data plot shows the measured wind speeds. Note that the 
dynamical flight mode leads to an air path speed of factor 3^ higher than 
the estimated wind speed at flight altitude which is significantly higher than 
the wind speed measured aboard the ship. 



we wanted a clear outline of the major features of the 
algorithm. They are now briefly summarized to give a more 
comprehensive understanding. 

First of all we have completely excluded the ship from our 
treatment by arguing that consideration of the apparent wind, 
which is the wind speed with respect to the ship's motional 
frame of reference, is an appropriate approximation for the 
basic design of the control system. Furthermore the choice of 
pattern parameters |25|, the optimization of the towing force 
with respect to the ship forward direction lIZTl as well as 
influences on dynamics due to waves are improvements to 
an operational kite autopilot but have not been presented here. 

It is worth mentioning that both, the thorough choice of 
the sensor set up and preprocessing algorithms, contribute a 
significant and crucial part to the presented controller perfor- 
mance and robustness. An important task is the estimation 
of the angle if) which involves not only sophisticated inertial 
navigation algorithms but also takes into account estimates 
for wind speed and wind direction at flight altitude as these 
may vary on the timescale of minutes. A further discussion of 
these mainly technical and cumbersome topics will be subject 
to future publications. 

In our discussions we have made use of constant and long 



line lengths which is the common operating point of our 
system. Winching of the towing line is done while launching 
and recovering our system and goes along with extra dynamic 
effects which are considered as additional correction terms 
to the presented equations. Especially for control during 
launch and recovery at shorter line lengths we use the idea 
of receding horizon feedback from the predictive control 
(NMPC) approaches f28l in order to provide ifj^. However 
our optimization computation is done analytically (thanks to 



the special and simple structure of the design model, see III I 
as opposed to having a numerical solver. 

Wind gusts may also go along with downwinds leading to 
a free flight situation of the usually constrained system. This 
situation involves a completely different dynamics. While (jT]) 
still holds in large part the kinematic changes significantly. It 
is a crucial advantage of the particular choice of coordinate 
system (see section [llji that controller input values behave 
in a favorable way during these exceptional occurrences. We 
would like to note that also disturbances due to excitations of 
internal modes of the real system are effectively suppressed 
although these modes are not considered explicitly in the 
'design model'. 

A more elaborate challenge is the avoidance of and response 
to stall situations. Systematic experimental tests on this are 
hardly feasible. Nevertheless it is an important issue and 
subject to current research and development activities. 

VIII. Summary 

In this paper we have presented a simple dynamical model 
for the dynamics of constrained kites. A key point has been 
that for controller design the complete aerodynamics can be 
reduced to a law involving only one or two quantities (see (j?]) 
and (|9]l, respectively). We have discussed flight data to justify 
that this reduction describes reality to a surprising high degree. 
Utilizing this we have developed a cascaded controller, based 
on the model following principle approach, and proved the 
effectiveness with flight test data. Finally the main principles 
of flight pattern generation have been given. 

We would like to finish this paper by emphasizing that the 
presented results hold for the SkySails towing kite system but 
the underlying equations are commonly valid for the dynamics 
of constrained kites. Thus a lot more applications of the 
presented controller are feasible especially in the fascinating 
upcoming field using kites in order to open up the green 
resource of wind power. 

Appendix 
Derivation of System Dynamics 



The system vectors read (see Fig. |2]i: 



eroll 



Spitch" 



— sin d cos ij) 

- cos Lp sin -0 + sin (p cos d cos ^ 
■ sin Lp sin tp — cos (p cos d cos ^ 

sin 1? sin 

- cos Lp cos — sin ip cos i? sin ■0 

- sin Lp cos "0 + cos ip cos d sin ^ 



(12) 



(13) 
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— COS \ 

— sin(/5sin-(? . (14) 
cos ip sin'd J 

As discussed in section|Il]we neglect gravity. Thus the problem 
becomes independent of (p and can be treated for = 
without loss of generality which implies the following basis 
vectors: 



€!pitch~ 



--yaw 



— sin d COS 

— sinip 

— cos 1? cos ip 

sin I? sin ip 

— cos V' 
cos "d sin ^ 

— cos ^ 


sini? 



(15) 



(16) 



(17) 



The air flow Va of the flying system is 




^roll ^pSpitch ■ 



(18) 



The first term describes the external wind, the subsequent two 
terms the flow due to the kinematic speeds Vr and Vp with 
respect to the basis vectors eioii and epitch- 

Considering the basic aerodynamics of an airfoil [22J as 
a very simple model we claim validity of the following two 
conditions: 

1) The airflow vector lies in the ej-oii-Leyaw -plane which 
means 

(epitch, S'a) -0. (19) 

2) The airflow direction with respect to e^ow and eyaw is 
given by the glide ratio E which is the ratio between 
lift and drag coefficients E = L/D, 



E. 



(20) 



Insertion of the definitions ([T5]l-([T7]i into ( 19 1 and ( [20) i yields 
the velocity components 



vo 



= sin 1} sin ?/> 

— iJcosi? — sin 1? cost/). 



(21) 
(22) 



The airflow in roll direction Va, measured by the anemometer 



in the control pod, can be calculated using ( 15 1 and ( 18 1 by: 

Va ^ ~ {va,eroii) ^ voEcos-d. (23) 

Geometric and kinematic considerations lead to the relation 
1 



'd = — {vr cos Ip — Vp sin tp) 



(24) 



Using (21 1, (22 1 and ( [23] l we get for the dynamics of d the 
following differential equation: 



d = 



COSTp - 



tan I? 
E 



(25) 



Similarly we obtain for ip the equation 
1 



Lsini? 



{—Vr sin Ip — Vp cos Ip) : 



and using (21 1, (22i and (23 i the equation of motion 



ip 



L sin?? 



sinip. 



(26) 



(27) 



For constant values of ^ one obtains the steady-state solution 
of ^ as 

'doW = arctan(£:cosV')- (28) 
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